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Abstract. Characteristic points have been a primary tool in the study of a gen- 
erating function defined by a single recursive equation. We investigate the proper 
way to adapt this tool when working with multi-equation recursive systems. 

1. Introduction and Preliminaries 

Recursively defined generating functions play a major role in combinatorial enu- 
meration, for example (see [7]) when studying lattice animals, context free languages, 
random walks on free groups, directed walks in the plane, colored trees, and Boolean 
expression trees. The subject started in 1857 with Cayley's [1] fundamental equation 
for the number t„ of rooted unlabeled trees: 



(1) = X 



n>l n>l 



He used this to recursively calculate (with some errors) the first dozen values of 
tn, and later applied his method to recursively enumerate certain kinds of chemical 
compounds. 

The generating function for trees is T{x) = J2n>i^nx"' . In 1937 Polya (see pTO] 
converted ([1]) into 

(2) T{x) = X ■ exp 

m>l 

a form to which he was able to apply analytic techniques to find asymptotics for the 
tn, namely he proved 

(3) tn ~ Cp-"n-=^/2 
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where p is the radius of convergence of T(x), and C a positive constant^ A similar re- 
sult held for the various classes of chemical compounds studied by Cayley. Although 
the function T{x) was not expressible in terms of well-known functions, nonetheless 
Polya showed how to determine C and p directly from ([2]). Polya's methods were 
applied to nearly regular classes of trees in 1948 by Otter [15]. 

In 1974 Bender [T], following Polya's ideas, formulated a general result for how to 
determine the radius of convergence p of a power series T{x) defined by a functional 
equation F(a:, y) = 0. Bender's hypotheses guaranteed that p was positive and finite, 
and that r := T(p) was also finite. His method was simply to find (p, r) among the 
solutions (a, b) (called characteristic points) of the characteristic system 



A decade later Canfield |3] found a gap in the hypotheses of Bender's formulation 
when there were several characteristic points. In the case of a polynomial functional 
equation, Canfield sketched a method (based on analyzing the singularities of an 
algebraic curve and how the branches relate to one another) to determine which of 
the characteristic points gives the radius of convergence of the solution y = T(x). 

In the late 1980s Meir and Moon [13J focused on a special case of Canfield's 
work, namely when F{x,y) = is of the form y = G{x,y), where G{x,y) is a 
power series with nonnegative coefficients. The interesting cases were such that 
setting T{x) = G{x, T{x)), with T{x) an indeterminate power series, gave a recursive 
determination of the coefficients of T{x). One advantage of their restricted form of 
recursive equation was that there could be at most one characteristic point. This 
formulation was adopted by Odlyzko in his 1995 survey paper [Tl] as well as in 
the recent book [7J of Flajolet and Sedgewick. These publications have focused on 
characteristic points in the interior of the domain of convergence of G{x,y), in the 
context of proving that p is a square root singularity of the solution y = T{x). If 
(p, r) is on the boundary of the domain of G{x, y) then p may not be a square-root 
singularity of T(x). 



F{x,y) 







dF 

dy 



{x,y) 



0. 



In ^ wc found this law so ubiquitous among naturally defined classes of trees defined by a 
single equation that we referred to it as the universal law for rooted trees. 
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Most of the areas of application listed in the first sentence of this introduction 
actually require a recursive system of equations 



written more briefly as y = G{x, y). (A precise definition of the systems considered 
in this paper is given in ^) This rich area of enumeration has been rather slow in it 
development. In the 1970s Berstel and Soittola (see [7] V.3) carried out a thorough 
analysis of enumerating the words in a regular language using recursive systems of 
equations that were linear in yi, . . . , y^- However it was not until the 1990s that pub- 
lications started appearing that used multi-equation non-linear systems. Following 
the trend with single recursion equations y = G{x, y), the focus has been on systems 
y = G(a;,y) where the (^^(Xjy) are power series with non-negative coefficients. 

In 1993 Lalley [10] considered polynomial systems in his study of random walks on 
free groups. In 1997 Woods [17] used one particular system to analyze the asymptotic 
densities of monadic second-order definable classes of trees in the class of all trees. In 
the same year Drmota [5] extended Lalley's results to power series systems. Lalley's 
and Drmota's results were for a wide range of irreducible systems, that is, systems 
in which each variable yi (eventually) depends on any variable yj. An irreducible 
system of the kind they studied behaves in some ways like a single equation system, 
for example, the standard solution yi = Ti{x) is such that all the Ti{x) have the 
same finite positive radius p, the Xj := Tj(p) are all finite, and the asymptotics for 
the coefficients of Tj(x) is of the Polya form Cip~^n~^^'^. 

Thus, as has been the case with single equation systems, it is desirable to find the 
radius of convergence p even though the solutions Ti{x) may be fairly intractable. 
The natural method was to extend the definition of the characteristic system from a 
single equation to a system of equations by adding the determinant of the Jacobian 
of the system, set equal to zero to, to the original system. The solutions of such a 
characteristic system will again be called characteristic points. 

Under suitable conditions one can find (p, r) among the characteristic points. To- 
date, however, the necessary study of characteristic points (a, b) for systems, so 
that one can locate (p, t), has been essentially non-existent. Filling this void is the 
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major goal of this paper. In December, 2007, one of the authors discovered that in 
the systems studied by Drmota it was possible that there could be more than one 
characteristic point — this was communicated to Flajolet and appears as an example 
in [7] (p. 484). The main objective of this paper is to give conditions to locate (p, t) 
among the characteristic points, if indeed (p, t) is a characteristic point. A review 
of and improvements to the theory of the single equation case (see Proposition [2^ 
are also given. It turns out that even if there is a characteristic point of a system 
y = G(x, y) in the interior of the domain of G(x, y), one cannot claim that the 
asymptotics for the coefficients of the solutions Tj(x) will be of the above Polya form 
see Examples [lEl [13)1 

Before proceeding with the study of such systems some background material is 
need. 

1.1. The extended nonnegative real numbers. 

Extend the usual operations on [0, oo) to [0, oo] in the obvious way as follows: 

c + oo = oo for c G [0, oo] 
c ■ oo = oo for c G (0, oo] 

the usual infinite sum if all c„ G [0, oo) 
oo if some c„ = oo. 



n 

Here the usual infinite sum is oo if the series diverges 



1.2. Formal power series in several variables. This section gives the essential 
definitions that lay the foundations for working with formal power series in several 
variables. The standard number systems are: 

^In 1997 Drmota fS^ appears to claim that having a characteristic point in the interior of the 
domain would lead to Polya asymptotics — however these examples show this not to be the case. 
In his 2009 book [B] this hypothesis is replaced with one regarding minimal characteristic points, 
which seems somewhat at odds with our Lemma [251 which says that the characteristic points form 
an antichain with the characteristic point (a, b) of interest having the largest value of a among 
the characteristic points. Theorem [36] of §6.11 is a restatement of Drmota's result which makes it 
clear which characteristic point is of interest, namely the one (if it exists) such that the Jacobian 
of G(x,y) has 1 as its largest real eigenvalue. 
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the set N = {0, 1, . . .} of nonnegative integers, the set Q of rational 
numbers, the set IR of real numbers, and the set C of complex numbers. 
For the hnearly ordered set R of real numbers one has the posets of real-valued 

functions on X, where the partial ordering is given by / < if f{x) < g{x) for all 

X & X. Familiar examples are: 

(a) n-vectors y — {vi, . . . , by setting X — {1, . . . ,ri\ 

(b) m X n-matrices M, by setting X — {1, . . . , m} x {1, . . . ,n} 

(c) formal power series in k-variables A{xi, . . . ,Xk) by setting X — N'^. In this 
case a function a from N'^ to R provides the coefficients, and one writes 

A(x) 5]a(i)x= 

A matrix (or vector) M of real numbers is non-negative (written M > 0) if each entry 
is non-negative, and positive (written M > 0) if each entry is positive. A power series 
A{'x) is non-negative (written >l(x) > 0) if each coefficient is non-negative. 

1.2.1. Composition of formal power series. 

For power series ^(u^i, . . . , Wm) and B({x), 1 < i < m, where the constant term 
of each Bi is zero, that is, 6^(0) = 0, define the formal composition 

C(x) ^(Si(x),...,S„(x)) 

by defining the coefficient function as follows: 

c(i) := 5^[x'] a{i) ■ B,{^y^ ■ ■ ■ Bmi^y- 
j>o 

Requiring that the constant term of the Bi{^) be guarantees that for each i only 
finitely many terms in this sum are nonzero. Consequently C(x) is indeed a formal 
power series. 

1.2.2. The function defined by a formal power series. 

A power series A(x) in k variables defines a partial function, also denoted vl(x), 
on M'^ (or C^) by setting 

(5) ^(c) := ^ ^ 

n>0 iiH hifc = n 
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whenever the sum converges. 

For A{x) a nonnegative power series in k variables, and for c G [0, oo]'^, A{c) = oo 
if the series ([5]) diverges, that is, if 



When one speciahzes to nonnegative power series, and focuses on nonnegative argu- 
ments c, then the general study of power series functions A{x) becomes much easier, 
in part because one does not have to worry about rearrangements in ([5]) affecting 
convergence, or the value to which it converges, but also because one is dealing with 
a function which is continuous, and is monotone nondecreasing in each variable. A 
nonnegative power series A{x.) in k variables defines a left-continuous function from 
[0, oo]'' to [0, oo], that is, if h~ in [0, oo]'' then A{aj) A(h). 

1.2.3. The derivatives of a formal power series. Derivatives of [nonnegative] formal 
power series give [nonnegative] formal power series: 



Working with the formal derivative of formal power series instead of the analytical 
derivative will impact the discussion of characteristic points. To see why this might be 
the case just note, for example, that the power series A{x) = Y^x'^/n^ has the formal 
derivative A'{x) = ^x^'^/n^ which converges at a; = 1, whereas the analytical 
derivative does not exist at x = 1. 

1.2.4. Holomorphic functions and a law of permanence. 

A complex- valued function /(x) of several complex variables is holomorphic at c 
if it is continuous and differentiable in a neighborhood of c. The notation [a, b] is 
short for 






j<n iiH hifc = j 




The notation A^^ is also used for the partial derivative dA/dxj. 



[ai,6i] X ■ ■ ■ X [ak,bk\. 



Proposition 1 (A Law of Permanence for Functional Equations). 

Suppose A(x), -B(x, y) > 0. If there is an e > such that 

A(x) = fi(x,A(x)) < oo /orxG[0,e] 
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then 

A(x) = fi(x,A(x)) /orxG [0,oo]. 
If furthermore a > and A(a) < oo then 

A(x) = 5(x, yl(x)) /or < Oi, l<i<fc 
and A(x) holomorphic for \xi\ < ai, 1 < i < k. 

Proof. This is a special case of Hille's law of permanence for functional equations 
given in §10.7 of Vol. 2, [9]. □ 

1.3. The Perron-Probenius theory of nonnegative matrices. The key to the 

main results of this paper are some simple observations based on the well-known 
Perron- Frobenius theory of nonnegative matrices that was developed ca. 1910. 

Proposition 2. Let M be a nonnegative nonzero k x k matrix with real entries. 

(a) M has a real eigenvalue. 

(b) The largest real eigenvalue A(M) is positive and is given by 

AIM) = max mm . 

x>0 l<i<fc Xi 

(c) A(M) is a simple root of the characteristic polynomial pm{X) = det(A/ — M). 

(d) The eigenspace belonging to A(M) is 1- dimensional, generated by a unique 
positive normalized eigenvector \m ■ (Normalized means the sum of the entries 
IS 1). 

Proof. (See §2 of Gantmacher [8].) □ 

Note that Proposition [2](b) implies that for some x > one has A(M) equal to 

(Mx), 
mmi<i<fc . 

X'l 

Corollary 3. 

(a) A positive k x k matrix M , k >2, has all diagonal entries < A(M). 

(b) A(X) is a nondecreasing function on the set of nonnegative matrices, that is. 
Ml < M2 implies A(Mi) < A(M2). Furthermore if every row [column] sum 
of Ml is less than the corresponding row [column] sum of M2 then A(Mi) < 
A(M2). 
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(c) A(X) is a continuous function on the set of nonnegative matrices, where the 
matrices are thought of as points in k"^- space. 

Proof. Since Mvm = ^{M)wm one has for 1 <i <k, 

< maVAiii) < ^ mijVMij) = A{M)vM{i), 

since each mijVMU) > 0- Thus ma < A(M). This proves (a). 

For (b) suppose < Mi < M2. Note that for x > and 1 < i < k one has 

(Mix), ^ (M2X), 



Xi Xi 

By Proposition [2](b) it follows that A(Mi) < A(M2). If in addition every row sum of 
Ml is less than the corresponding row sum of M2 then for x > and 1 < i < k 

(Mix)i ^ (M2X), 



Xi Xi 

With this inequality Proposition [2](b) shows that A(Mi) < A(M2). 

Since a square matrix M and its transpose M* have the same characteristic poly- 
nomial, and hence the same eigenvalues, one can replace the condition on row sums 
by one on column sums, finishing the proof of (b). 

To see that (c) is true note that the largest real eigenvalue of M is simple, so the 
graph of the characteristic polynomial PmW crosses the axis at A = A(M) with a 
non-zero slope. Consequently this crossing point changes continuously with changes 
in the entries of M. (Note: A special case of item (c) is stated on p. 2103 of Lalley 
[TU] . for certain Jacobian matrices denoted evaluated along certain curves.) □ 

2. Well-conditioned systems 

The next definition gives a version of essentially well-known conditions which en- 
sure that a system y = G{x,y) has power series solutions i/i = Ti{x) of the type 
encountered in generating functions for classes of trees. (See Drmota [5], [6].) 

Definition 4. A system y = G(x, y) is well-conditioned if it satisfies 

(a) each Gi{x,-y) is a power series with nonnegative coefficients 

(b) G{x, y) is holomorphic in a neighborhood of the origin 

(c) G(0,y) = O 
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(d) for alii, Gi{x,0) 

(e) det (/ — Jg(0, 0)) 7^ where Jg is the Jacobian matrix I — 

(f) the system is irreducibl^ 

d'^GAx, y) 

(g) for some i,j, k, — „ ' ' — ^ (so the system is nonlinear in y). 

dyjdyk 

Remark 5. One can replace condition (b) by (b'); G(x, y) converges at some positive 
(a,b). 

2.1. Solutions of Well-Conditioned Systems. 

Proposition 6. Ify = G(x, y) is a well- conditioned system then the following hold: 

(i) There is a unique vector T(x) of formal power series Ti{x) with nonnegative 
coefficients such that one has the formal identity 

(6) T(x) = G(x,T(x)). 

(ii) Equation ([6]) gives a recursive procedure to find the coefficients of the Ti{x). 
(ill) Equation ([6]) holds for x G [0,cx3]. 

(iv) All Ti{x) have the same radius of convergence p G (0, 00) and all Ti{x) con- 
verge at p, that is, Xj := Tj(p) < 00. 

(v) Each Ti{x) has a singularity at x = p. 

(vi) If {p,t) is in the interior of the domain of G{x,y) then 

det {I - Jg{p,t)) = 0. 

Proof. By conditions (b) and (e) of Definition HJ tlie implicit function theorem guar- 
antees a unique solution y = T(x) to the system in a neighborhood of the origin. 
This solution is holomorphic in a neighborhood of the origin, so each Tj(x) has a 
radius of convergence > 0. Let p be the mapping on m-tuples of nonnegative power 
series defined by p[A(x,y)) = G(x, A(x, y)). In view of Conditions (c) and (d) of 
Definition m p'^{x, 0) gives a sequence of m-tuples of nonnegative power series that 
converges to an m-tuple of nonnegative power series that provides a formal solution 
to the system, so it must be T(a;). Then by Proposition [1] the equation 

T(x) = G(x,T(x)) 



This means the non-negative matrix Jg is irreducible. 
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holds for X G [0, oo]. Since the system is irreducible, if x > is such that some 
Tj(x) = oo then all Tj(x) = oo. This means all Tj(x) have the same radius of 
convergence p which, as noted above, must be > 0. From the nonlinearity condition 
(g) of Definition m it follows that p < oo and all = Tj(p) < oo. Each Ti{x) has a 
singularity at x = p by Pringsheim's Theorem, which gives (v). Thus if (p,t) is in 
the interior of the domain of G then the Implicit Function Theorem guarantees that 
the determinant of the Jacobian of the system must vanish, giving (vi). 

□ 

The sequence T(a;) of power series described in Proposition [6] is the standard 
solution of the system, and the point (p,t) is the extreme point (of the standard 
solution, or of the system). From (jS]) one has T(0) = 0, so the standard solution 
goes through the origin. The set 

Dom+(G) := {(a, b) : a, 6i, . . . , 6m > and Gj(a, b) < oo, 1 < i < m} 

is the positive domain of G. For (a, b) G Dom+(G) let 

A(a,b) := A(jG(a,b)), 

the largest real eigenvalue of the Jacobian matrix Jc,{a,h). 

2.2. Characteristic Systems, Characteristic Points. Flajolet and Sedgewick 
[7] VII. 6 define the characteristic system of (jlj) to be 

yi = Gi{x,yi,...,ym) 

ym = Gm{x,yi, . . . ,ym) 

^ = det(/- JG(x,y)). 

Let the positive solutions (a, b) G to this system be called the characteristic 

points of the system]^ Requiring that (p, t) be a characteristic point in the interior 
of the domain of G(x, y) has been crucial to proofs that x = p is a square-root 
singularity of the Tj(x), leading to the asymptotics tj(n) ~ Cjp~"'n~^/^ for the 
non-zero coefficients. 

'^Flajolet and Sedgewick ([J Chapter VII p. 468) only consider characteristic points in the interior 
of Dom+(G). 



CHARACTERISTIC POINTS 



11 



There is considerable interest in finding practical computational means of estimat- 
ing p — for then, at least in the case that (p,t) is in the interior of the domain of 
G, one has a good estimate for the growth of the non-zero coefficients of the Ti(x), 
namely up to a constant it is p~^''n~^/'^. (To refine this further one needs estimates 
of the constants Cj. ) 

For the case that the Gj(x, y) are polynomials we know that (p,t) will be among 
the characteristic points and in the interior of the domain of G. However until 
now, even in the polynomial case, no general attempt has been made to characterize 
(p, r) among the characteristic points of the system — with one exception, namely the 
1-equation systems. 

3. A Collection of Basic Examples 

The following examples explore the behavior of characteristic points of well-conditioned 
systems. 

3.1. 1-equation systems. For a 1-equation system y = G{x,y) there can be at 
most one characteristic point, and if such exists it must be (p, r) — see Proposition 
[291 Thus if there is a characteristic point then one has, in theory, a way to determine 
(p, r), namely find the positive solution to the characteristic system. If there is 
no characteristic point then, at present, there is no simple alternative method to 
determine (p, r). 

In general (p, r) lies either in the interior of the domain of G or on the boundary 
of this domain — in either case x = p is a singularity of the standard solution T{x). 
If (p, r) is in the interior of the domain of G then x = p is a square-root singularity 
of T{x). The possibilities for the nature of this singularity when (p, r) is on the 
boundary of the domain of G have not been fully classified. 

3.2. Examples for 1-equation systems. For 1-equation systems the following two 
examples show the three kinds of behavior described above, namely: (i) there is a 
characteristic point which is an interior point and thus equal to (p, t), (ii) there is 
a characteristic point which is a boundary point and thus equal to (p, t), and (iii) 
there is no characteristic point. 

Each example starts with an equation y = G{x, y) where the characteristic point 
(p, r) is in the interior of the domain of G{x,y). Then the example is modified to 
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give a system y = G*{x, y) with (p*, r*) on the boundary of the domain of y). 
{p*,T*) is a characteristic point in Example [7] but not in Example [8l 

Example 7. Let G{x,y) = x{l For the characteristic system 

( y = + 1/2) 
\ 1 = 2xy 

of y = G{x,y) one has the characteristic point (1/2,1), an interior point of the 
domain of G{x,y), so for the standard solution y = S{x) of y = G{x,y) one has 
{p,t) = (1/2,1) . The established theory for such a system (see Chapter VII) 
shows that S{x) has a square-root singularity at x = p. 

Next let G*{x,y) = S{x){l + y'^)/2. For the characteristic system 

(y = S{x){l + y')/2 
I 1 = S{x)y 

once again the characteristic point is (1/2, 1), but now it is a boundary point of the 
domain ofG*{x,y). An examination of the standard solution (see Proposition [77]) 
of y = G*{x,y), namely y = T{x) = S(^S{x)/2), shows that it has a fourth-root 
singularity at x = 1/2. 

Example 8. Let G{x, y) = x{l + 2y + 2?/^) . The characteristic system 

( y = x{l + 2y + 2y^) 
I 1 = 2x{l + 2y) 



of y = G{x,y) has the characteristic point 

V2-1 V2 



an interior point of the domain of G{x,y), so for the standard solution y = S{x) 
of y = G{x,y) one has p = [V2 — l)/2 and r = V2/2. S{x) has a square-root 
singularity at x = p. 
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Next let G*{x, y) = a;(l + S{x) + y + 2?/^) . The standard solution of y = G*{x, y) 
is again y = S{x), so {p*,T*) = {p,t). The characteristic system 



( y = x{l + S{x) + y + 2y^) 
I 1 = x(l + 4?/) 



of y = G*{x,y) has no characteristic point since the only candidate is {p,t) and 



(p, r) is a boundary point of the domain of G*{x, y) whose location is not detected by 
the method of characteristic points. 

Remark 9. On p. 83 of their 1989 paper [13j Meir and Moon offer an interesting 
example of a 1-equation system without a characteristic point, namely y = y4(x)e^ 
where A{x) = (1/6) Xln^"/''^^- -^^^ characteristic system is 



so a characteristic point {a,b) must have 6=1, A{a) = 1/e. But 1/e is not in the 
range of A{x), so there is no characteristic point. One can nonetheless easily find 
{p,t) in this case since {p,t) must lie on the boundary of the domain of A{x)ey . 
Thus p = l, and then r = A(l)e^ = (7rV36)e^, so t ^ 0.41529. 

The paper goes on to claim that by differential equation methods one can show that 
the standard solution y = S{x) has coefficient asymptotics s{n) G/n. However this 
cannot be true since such a solution would diverge at its radius of convergence p = 1 
(see [2\), whereas the given equation y = A{x)e'^ is nonlinear in y, so the solution 
must converge at p. 

3.3. 1-equation framework. This subsection gives a framework for 1-equation ex- 
amples which will be useful for building the 2-equation examples in §3.41 

Proposition 10. Let A{x) be the standard solution of the recursive quadratic equa- 
tion 



p{l + At) = (1/2) {V2 - 1) (1 + 2V2) ^ 1. 



y = A{x)ey, 1 



v4(x)e^ 





where a > and b > 0. Then the following hold: 
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(a) 



(9) 



A{x) 



1 



(1 - ax) - - axy - 4:bx^ 



) 



2bx 



(b) A{x) has non-negative coefficients. 

(c) A sufficient condition for A{x) to have integer coefficients is that a and b are 
integers. 

(d) A{x) has a positive radius of convergence pA given by 

0) PA = 

a + 2Vb 

(e) ta '■= A{pa) is finite and is given by 



(f) Pa is a square-root branch point of the algebraic curve defined by (IHD- 

(g) {pa, Ta) is the unique characteristic point of that is, it is the unique positive 
solution {x, y) to 



Proof. Writing ([8]) in the form bxA^ + (ax — l)y + x = 0, there are two solutions: 



The standard solution goes through the origin, which is ([9]). (The other solution 
blows up at the origin.) 

(b) and (c) are standard facts about recursively defined functions. 

The radius of convergence pA is found by setting the discriminant equal to 0, that 



1 



y = x(l + ay+ by^) 
1 = x(a + 2b?/). 





Then from and (ITU]) one calculates ta := A{pa)- 

^ 1 - apA ^ J_ 
~ 2bpA ~ Vb' 
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Equation ([8]) is a polynomial equation, so it defines an algebraic curve A*{x). For 
such curves the only singularities are poles and branch points, and pa is clearly of 
the latter kind. 

(g) is a straightforward calculation. □ 

Let S{x) be the standard solution to the quadratic equation 
(11) y = x(l + C2/ + V) 

where c > and c) > 0. Then, from the previous result, 

S{x) = ^ ((1 - ca;) - V(l - cx)2 - 4dx^^ 

(P5, Ts] ^ 



{ps,Ts) is the unique positive solution {x,y) to 

y = xil + cy + dy"^) 
(12) 1 = x{c + 2dy), 

and ps is a square-root branch-point singularity of the algebraic curve defined by 

m. 



Proposition 11. Given a, c > and b,d > let A{x) be the standard solution of 

y = x{l + ay+by"^) 
and let S{x) be the standard solution of 

y = x{l + cy + dy'^). 
Let T{x) be the standard solution of the recursive quadratic equation 
(13) y = A(x)(l + C2/ + V). 

Then the following hold: 

(a) T{x) = S{A{x)). 

(b) Tix) = ((1 - cAix)) - v/(l - cA{xW - 4M(a;)2) . 

(c) T{x) has non-negative coefficients. 

(d) A sufficient condition for T{x) to have integer coefficients is that a, b, c,d are 
integers. 



16 JASON BELL, STANLEY BURRIS, AND KAREN YEATS 

(e) // = c + 2v^ then 

and T{x) has a fourth-root singularity at px- 

Proof. Since S{x) solves ( fTTll . 

S{x) = x{l + cS{x) + dS{xy). 

Substituting A{x) for x gives 

SiA{x)) = A{x){l + cS{A{x))+-OS{A{x)f). 

Thus T(x) := S{A{x)) is the standard solution to y = A{x){l + cy + dy"^). 

Since ([8]) defines an algebraic curve with a square-root branching point at pa, and 
likewise (fTTj) defines an algebraic curve with a square-root branching point at ps, the 
composition T{x) = S{A{x)) defines part of an algebraic curve with a fourth-root 
branching point Pa- D 

Given the premises of Proposition [HI the restriction ^/b = c + 2y/d is called the 
critical composition condition (CCC); this is the condition needed for T{x) = 
S{A{x)) to be a critical composition (as defined by Flajolet and Sedgewick p. 411). 

3.4. Multi-equation systems. For systems y = G{x, y) of two or more equations 
one again has (p,t) in the interior or on the boundary of the domain of G. However 
in general the situation is more involved than in the case of a single equation. Unlike 
the one-equation case, it is possible that there are multiple characteristic points, in 
which case they form an antichain (see Lemma 1261) . If (p,t) is a characteristic point 
then (see Proposition [31!) among the characteristic points (a, b) it is the one with 
the largest value of a. 

The goal of this paper is to show that there is a simple additional constraint that 
one can put on characteristic points so that they are just as useful in the multi- 
equation case as in the single equation case. Namely (see ^ the property 

1 is the largest positive eigenvalue of JG{a.,h) 

gives a definitive test for a characteristic point (a, b) to be {p,t). A characteristic 
point with this property will be called an eigenpoint. 
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Before considering eigenpoints, a family of 2-equation systems is constructed, based 
on the results of §3.3[ These systems yield the main examples in this paper, namely 
Examples [I6] and [T71 

Proposition 12. Suppose 

a, ci>0, b,C2,c)>0, Vb = c + 2\/d, c = Ci + C2. 

Let A{x) be the standard solution of 

y = x{l + ay+by^), 
let S{x) be the standard solution of 

y = x{l + cy + dy^), 
and let T{x) be the standard solution of 

y = A{x){l + cy + -Qy^). 

Then the following hold: 
(a) The quadratic system 

{SYS) : 



A{x){l + ciT{x) + C2y2 + 'Qyi^) 

A(x)(l + CiT(x) + C2?/l + C)l/2^) 



is well- conditioned, and the standard solution is yi = y2 = T{x) 
(b) The extreme point (p, r, r) of (SYS) is given by 



It is on the boundary of the domain of (SYS). 

(c) T{x) = S{A{x)) has a fourth-root singularity at x = p. 

(d) A positive point {x, y, y) is a characteristic point of {SYS) iff either 

A{x) (^C2 + 2^c)(l + CiT(x)) ^ 

1 - C2A{x) 

2dA{x) 



(*) 
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or 



Aix) (^C2 + 2y/c22 + 0(l + CiT(x)) ^ 

1 + C2A{x) 



' 2dA{x) ■ 

(e) // Ci = then there are exactly two characteristic points of the form {x,y,y): 
the first is {p,t,t), a boundary characteristic point obtained from and 
the second is the unique positive solution to an interior characteristic 
point. This is the only case where (^) contributes a characteristic point, namely 
{p,t,t), and this is the only case where {p,t,t) is a characteristic point. 

(f) //O < Ci = 2c2 then there is a unique characteristic point of the form {x,y,y): 
it is the unique positive solution to (irk) and it is a boundary point different 
from (p, r, r) . 

(g) //O < Ci < 2c2 then there is a unique characteristic point of the form {x,y,y): 
it is the unique positive solution to (irk) and it is an interior point that is 
different from (p, r, r) . 

(h) //2c2 < Ci then there are no characteristic points of the form {x,y,y), so again 
(p, r, r) is not a characteristic point. 

(i) The second characteristic point in (e) and the unique characteristic points in 
(f) and (g) are given explicitly by 



X 



where 



ac + 2c2 + f + b + (a + 2c)yFTf 
C + C2 + Vc^ + f 
2c) 



f = -6ciC2 + 3C2 + 4c). 



Proof. The system is irreducible since C2 > 0. T{x) satisfies ( fT3i) . consequently 
Hi = y2 = T{x) is a solution of {SYS). This must be the standard solution of 
{SYS), so the extreme point of {SYS) is the stated (p, r, r) (see §3.3p . a point which 
is not in the interior of the domain of the system {SYS) since A{x) has a singularity 
at Pa, and p^ = p by (CCC). This takes care of items (a)-(c). 
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The characteristic system of (SYS) is 

f 2/1 = A(x)(l + CiT(x) + C22/2 + f2/i') 
(CSYS): I y, = A(x)(l + CiT(x) + 022/1 + 02/2^) 

[o = (l-20A(x)2/i)(l-20A(x)2/2) -C2M( 

so a positive point {x,y,y) is a characteristic point of S iff (x,?/) satisfies 



{CSYS') : 



y = A(a;)(l + CiT(x) + C22/ + V) 
= {l-2dA{x)yY -C2'A{ 



xf'. 



This gives the conditions (*) and (*^) of item (d) as follows. 
{CSYS') yields 

(14) 1 = {2Dy±C2)A{x), 

and then 

(2c)2/±C2)2/ = 1 + CiT(x) + C22/ + V, 



so 

.2 



(15) V+(-l±l)c22/ = l + CiT(x). 

From (fT^ and (fTSll there are two cases to consider: 
Case 1: 

1 = (2D2/ + C2)A(x) and Oy^ = 1 + CiT(a;). 

Case 2: 

1 = {2l3y - C2) A{x) and dy^ - 2c2y = 1 + Cir(a;). 

In each case, solving for y in the two equations and equating the results gives an 
equation to solve for x: 
For Case 1: 



(16) y 1-^2A(X) _ ,/l + Cir(x) 



2dA{x) 
For Case 2: 



(17) 
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Note that one cannot use the negative sign before the square root when solving for 
y in Case 2 since that would make y negative. 
CASE 1: 

To find candidates for x in the first case one needs to solve f|T6|) . which implies: 



1 = A{x)(^C2 + 2^d{l + CiT{x))^ 



By the monotonicity of the right side, there is at most one positive x solving this 
equation. If such exists, then one finds y by (fT6|) : 



Thus if (fTSj) has a positive solution for x then defining y by (fT9l) gives the positive 
point ix,y,y). 

Conversely, given a positive point (x, y, y) satisfying ( fT8|) and ( fT9l) . the above steps 
are reversible to show that (x, satisfies {CSYS'), so {x,y,y) is a characteristic 
point of (SYS). 

CASE 2: 

To find candidates for x in the second case one needs to solve ( JT7j) . which implies: 



(20) 



1 = A{x) (^C2 + 2^02^ + (1 + CiT{x)) ^ 



By the monotonicity of the right side, there is at most one positive x for this case. 
If such exists, then let 

Thus {x,y,y) is a positive point. 

Conversely, given a positive point (x, y, y) satisfying (120|1 and (12T!) . the above steps 
are reversible to show that (x, y) satisfies {CSYS'), so {x,y,y) is a characteristic 
point of (SFS). This finishes the proof of (d). 

For items (e)-(h) note that (p, r, r) is a characteristic point of (S'FS') iff the second 
equation of {CSYS') holds at this point, that is, iff 

(22) = {l-2dA{p)rY -C2^A{p)'. 
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Since t = ts and A{p) = ta = Ps, (1221) holds iff 

(23) = {l-2dpsrsY -C2''pI. 
Now from we have 1 = ps{c + 2dTs), so ([23]) holds iff 

= cVI-CaVI, 

and this holds iff = — iff Ci = 0. 
Let 

:= A(x)(^C2 + 2y/o(l + CiT(x)) ^ 

F2(a;) := (^C2 + 2^02^ + c) (l + CiT(a;)) ^. 

Then the first equations from (*) and (*T*r) are 

(24) 1 = Fi(x) 

(25) 1 = F^ix). 

The functions Fi{x) and -^2(0;) have radius of convergence = p, they are strictly 
increasing on [0,p], and -Fi(O) = -^2(0) = 0. Thus has a positive solution for a; 
iff 1 < Fi{p), and fl25l) has a positive solution iff 1 < F2{p). 

There are two possibilities for a positive solution x of 1 = Fi{x), for i = 1,2: 

(26) X = p is a solution iff 1 = Fi{p) 

(27) some x < p is a solution iff 1 < Fi{p). 



Now 



so 



Fi{p) = A(p) (^C2 + 2^0(1 + CiT(p))^ 

l = Fi(p) iff 1 = -_l-^(^c2 + 2y^7^^ 
iff ci = 0. 
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Likewise 

1 < iff ci^ < 0, 

which is impossible. 

Also, 



l = F2(p)iffl = ^^2v^ ' C2 + 2Vc2' + c) + CiVD 
iff Ci = 2c2, 



and 



1 < F2(p) iff Ci < 2C2. 

From the above one has the following cases, noting that Fi{x) < F2{x) on (0,p]: 

• ci = iff 1 = Flip) < Hp) 

• < ci = 2c2 iff Fi(p) < 1 =F2(p) 

• < ci < 2c2 iff Fi(p) < 1 < F2(p) 

• < 2c2 < Ci iff Flip) < Fiip) < 1- 

Finally now consider (i). With the assumptions of (e) we have Ci = 0. Condition 
(^) gives the characteristic point 



Consider The first equation of i-k-k) is independent of giving 

(28) Aix) = \ 

Solving ([9]) for x, using ( l28l) . and solving the second equation of i-k-k) for y using ( l28ll 
for y and rearranging gives the desired expressions for x and y in the special case 
ci = 0. 

Assume < Ci < 2c2. Substitute yi = T(x), y2 = Fix) into iSYS). Both 
equations of iSYS) give the same equation in Aix) and T(x), namely: 

(29) = -T(x) + A(x)(^l + (ci + C2)T(x) + c)T(x)'). 
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From the first equation of {-k-k) one also has 



(30) 



1 = A{x) + 2^cl + d(l + CiT{x)'^^ 



Solving (l30i) for T{x) gives a rational expression in A{x) which has denominator a 
multiple of A{x)'^ and has numerator a quadratic polynomial in A{x). Substitute 
this in for T{x) in ( !29l) . Simplifying this gives a quartic polynomial in A{x) equal to 
0. This polynomial can be solved explicitly for A{x). Let 



f = _6ciC2 + 3c2 + 4c) 



= 2ciC2 + 3c^ + 40 
If f 7^ then solving the polynomial gives 

-c + v/(ci + 2c2)2 + 40 _ -c + v/c^Ts 



A, 

A2 

A3 

A, 



2ciC2 + 3c^ + 4c) Q 



-c- V(ci + 2c2)' + 4c) ^ -c- y^Tg 
2ciC2 + 3c^ + 40 ~ 

-c + v/(ci - 2C2)^ + 40 ^ -c + y/^^ 
-6ciC2 + 3c^ + 40 ~ f 

-c- V(ci -2c2)2 + 40 -c-v/FTT 



-6C1C2 + 3ci + 40 f 
If f = then the quartic polynomial simplifies to the cubic polynomial 

1 - AcA{x) + A{cl + cl)A{xf + 16cciC2A(x)^ 

Solving for A{x) gives Ai and A2 as before and the third root is 

^ 2c 

All that remains to do is to identify which A^ gives the characteristic point and 
then solve for x from ([9]) and for y from the second equation of (**) and rearrange 
to obtain the explicit formulae. 

The identification of the correct Ai is given in the following Lemmas. 

□ 

In order to compare the Ai the following elementary lemma will be handy. 
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Lemma 13. Suppose c > 0, — <t<u,u^Q,t^Q. Then 



-c + ^JlF+l ^ -c + Vc^ + u 
t ~ u 

Proof. Simply calculate keeping track of the signs of t and u: 



-c + + t -> + Vc^ + M 



t u 



> , , if t, M < or t, u > 



<J lift<0<M 

o o o I > /7 T if t, < or t, M > 

M^t + t\ + 2cHu I ~ > 2tuJ{c^ + t){c^ + u) { 

<\ if t < < M 



^ M + t + 2c^ > 2v/(c2 + t)(c2 + m) 
^ (m - > 

The second and fourth implications are by squaring both sides and rearranging; in 
both cases both sides of the previous step are positive regardless of the signs of t and 
u. □ 

Lemma 14. With notation as in the proof of Proposition [iH z/f = then gives 
the characteristic point while if f ^ then A-^ gives the characteristic point. 

Proof. First determine when each Ai is positive. A2 is negative, and Ai and A^ are 
positive since Ci > 0, C2 > 0, and > 0. 

A^ is no smaller than Ai and hence positive by Lemma [T31 A4 has negative 
numerator and so is positive provided f < 0. 

Comparing the Ai 

< < < A3 (< A4 when f < 0) 

Also, can be discarded even when f < as it is greater than the radius, which 
due to (CCC) is l/(c + 2^^). This is because f + > so 



f+c2 + 2cV^ > > -(c + 2v^)v/c2 + f 
1 ^ -c-V^^f 



+ 2v^ 
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when f < 0. 

Another source of spurious solutions is the fact that (!29l) is quadratic in T. Solving 
( l29l) for T we are interested in the lower branch — the negative sign for the square 
root. So we must have 

1 - cA, - v/(l - cAjy - 4c)i^ _ 1 - 2c2Ai - (Scj + 4d)Aj 
2dAi ~ 4ciOA2 

Thus 

1 - cA,)2 - 4c)A^ = -1 + 2cA + {-2cl -jciC2 + Scj + 4D)Af 
Assuming Aj > leads to the requirement 

(31) -l + 2cAi + (-2c2-2ciC2 + 3c2 + 4c))A2 > 

As a polynomial in Aj this is negative at Aj = and has roots 



^ _ -Ci± v/(2C2-Ci)(2C2 + Ci)+4c) ^ -c± y/FT^ 

-2c2-2ciC2 + 3c2 + 4c) t) 

where 

f) = -2ci - 2ciC2 + 3cl + 4d 

When = (!3T!) becomes 

Ai > — 
- 2c 

Suppose f) = 0. Then f = 2ci(ci — 2c2) < and q = 2ci(ci + 2c2). Similarly to 
Lemma [T3l for any — < t, t 7^ 0, 



(32) 

since 

(132]) holds 




if t < 



■<=^ > by squaring both sides 

Applying the above with f and g we see that Ai does not satisfy the requirement 
but A3 does. Also, when f = 0, A5 trivially satisfies the requirement and so by the 
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same calculation Ai < A^. Thus A^ gives the characteristic point if f = and ^3 
otherwise. 

Now suppose P) 7^ 0. Note that since Ci < 2c2 we have f) > f, with equality when 
Ci = 2c2- Thus by Lemma [T3l 

A2<0<Ai< Crit+ < A3 (< when f < 0) 

When f) > 0, Crit_ < and so the positive part of ( !3TI) is below Crit_ and above 
Crit+. Also by fl32l) Crit+ < ^5, so if f = then A^ must give the characteristic 
point. On the other hand if f 7^ then A4 is out of the picture as it is either negative 
or larger than the radius and so A^ must give the characteristic point. 

Suppose f) < 0. Then f < 0. Similarly to Lemma [13] 

Crit_ > 

since 

-C- y/FT^ ^ -C- y/FTT 
^ - f 

^ (f)-f)c > -f)Vc' + f + fVc' + f)- 

The left hand side of the last line is nonnegative, so either the right hand is negative 
and we are done or we can square both sides and continue the chain of equivalences: 

which is true since the left hand side is nonnegative and the right hand side is 
negative. 

Thus in this case 

< < Ai < Crit+ < A3 < < Crit_ 

and the region where fl3T]) is positive is the region between Crit+ and Crit_. is 
larger than the radius and so again A3 must give the critical point. 

□ 

3.5. Examples of 2-equation systems. Now we look at three well-conditioned 
examples that show some of the varied behavior of characteristic points when one 
has more than one equation in the system. In the first example there are two char- 
acteristic points, both in the interior of the domain of G{x,y) and one of them is 
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(p, r). In the second example one has a characteristic point in the interior of the 
domain of G{x, y) and (p, r) is a characteristic point on the boundary of the domain. 
In the third example one has a characteristic point in the interior of the domain of 
G{x,y) but (p, t) is not a characteristic point. In the second and third examples, p 
is not a square-root singularity of the solutions. Such examples show the need for 
a more subtle use of characteristic points in the pursuit of information on (p, r) for 
multi-equation systems. 

Example 15. For the system of two equations 

Vi = x-{l+y2 + 2yl) 
y2 = X- {l + yi + 2yl) 

add 

{I - Axy^){l - Axy2) - x" = 

to obtain the characteristic system. This is a polynomial system, so all characteristic 
points will he in the interior of the domain; and since (p, ri, t^) is also in the interior 
it must he a characteristic point. Let (a, 6, c) he a characteristic point. Ifb^c then 
subtracting the second equation from the first gives 

b-c ^ a{c-b + 2b'^ -2c^). 

Dividing through by a{b — c) gives 

- = 2(6 + c)-l. 

a 

The third equation says 

substituting in the expression for 1/a from the previous equation gives 

(2c - 26 - 1) {2b - 2c - 1) = 1. 

Thus 

l-4(6-c)2 = 1, 

which is impossible ifb^c. Thus b — c, so the characteristic points are the positive 
triples (a, 6, b) satisfying 

b = a(l + 6 + 262) 
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= {l-4abf. 
From this the system has two characteristic points: 




^ (0.2612,0.7071,0.7071) 



2^3- 1 1 + ^3 1 + ^3 



(0.2240,1.3660,1.3660). 



Now we are left with determining which of the two characteristic points gives (p, ri, T2). 
By applying either Proposition [3l\ or LemmalSE below, it is the first of these. 

Next we apply the results of the previous section to create two examples showing 
that the behavior of the extreme point is relatively independent of whether or not 
there are characteristic points in the interior of the domain of G. 

Example 16. Let = 0, b = 9, Ci = 0, C2 = 1, and = 1. These numbers satisfy 
(CCC). Following the hypotheses of Theorem UB, let A{x) be the standard solution 
to y = x{l + 9y^) and consider the system 



Since ci = there are two characteristic points of the form (a, b, b) . The first is the 
extreme point 



which lies on the boundary of the domain, and the second is the interior point obtained 
from the formulas in Proposition ITM (i): 



Example 17. Let a = 0, b = 16, Ci = 1, C2 = 1, and = 1. These numbers satisfy 
(CCC). Following the hypotheses of Theorem UB, let A{x) be the standard solution 
to y = x(l + 16y^), and let T(x) be the standard solution to y = A{x){l + 2y + y"^). 
Consider the system 



yi = A{x) ■ {I + y2 + yl) 
7/2 = A{x)- (l + yi + yl). 



ip,n,r2) 



(1/6,1,1) 




yi = A{x)- {l + T{x) + y2 + yl) 
y2 = A{x)- {l + T{x)+yi+yl). 
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Since < Ci < 2c2, the extreme point 



(p,n,r2) 



(1/8,1 



1) 



is not a characteristic point, hut there is a characteristic point of the form (a, 6, h) 
in the interior of the domain of G given by the formulas of Proposition [TR - 



4. Characteristic Points of Well-Conditioned Systems 

From now on it is assumed, unless stated otherwise, that we are working with a 
well-conditioned system S : y = G(x, y) of m equations. 

4.1. Making substitutions in an irreducible system. A careful analysis of the 
characteristic points of S is easier if Jg{x, y) is a positive matrix for positive points 
(a;, y); this is the case precisely when no of entry JG{x,y) is 0. Fortunately there 
is a simple substitution procedure to transform the original system S into a well- 
conditioned system S* with 

(i) exactly the same positive solutions (a, b), and 

(ii) exactly the same set CV of characteristic points, 

even though the Jacobian Jq* {x, y) of the new system can be quite different from 
that of the original. 

The simplest substitutions to consider are ra-fold iterations G''"^ of the transfor- 
mation G defined by: 



This uniform substitution method is used in [7j (see p. 492) on aperiodic^ polynomial 
systems S. The following example shows that, in general (without the aperiodic 
hypothesis), this uniform method of substitution fails to produce a system S* such 
that the Jacobian matrix Jq* {x, y ) has all entries non-zero. However a more general 
form of substitution will provide the desired S*. 

well-conditioned system y — G{x, y) is aperiodic if the coefficients of each Ti(x) are eventually 
positive, T(x) being the standard solution — see |7|, P- 489. 



(a, b, b) 




G«(x,y) := G(x,y) 
G("+i)(a;,y) := G(")(a;,Gi(x,y),...,G„(a;,y)). 
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Example 18. Consider the irreducible system y = G(x, y) of 4 equations: 



yi 




■ ,yi) 


:= x(l+y2' + 2/4') 




= G2(x,yi,. 


■ ,yi) 


:= ^(l+yi^ + ya') 


2/3 


= Gz{x,yi,. 


■ ,yi) 


:= x{l+y^) 




= Gi{x,yi,. 


■ ,2/4) 





Lei M = Jg(") ■ r/ien is easy to check that Mn ^ iff' n is odd, and M12 iff' 
n is even. Thus for n > 1, JGi^ii^yj) has entries which are 0. 

One can transform S into a system S* where the Jacobian of G* has all entries 
non-zero by doing selective substitutions. For example, in the first equation of S 
replace one of the two y2 's by G2{x, y), giving the system 

yi = x{l+y2G2{x,y)+y4^) 
y2 = ^(l+yi^ + ya^) 
2/3 = x{l+y4^) 

^ yi = x{l+yi^) 

The first equation in this system is such that the right hand side depends on all 4 of 
the yi. Continuing in this manner one obtains a system 

{yi = G\{x,yi,...,yi) 
; 
7/4 = Q(a;,yi,...,7/4) 

in which every y) depends on each ofyi,..., ^4. 

Clearly every positive solution {a, b) of the original system E is also a solution of 
E*. We need to look at the question of whether or not there are any further positive 
solutions ofT,*, and how the set of characteristic points ofT,* compares to that ofE. 
One finds in general, with a little caution in choosing the substitutions (in order to 
maintain irreducibility), that one can create a S* that has exactly the same positive 
solutions, and the same characteristic points, as E. The need for a carefully selected 
sequence of substitutions yj ^ Gj{x,y) is largely removed in the following by using 
fractional substitutions yj ^ aGj{x,y) + (1 — o:)yj, a G (0, 1). 
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Given a system E : y = G(x, y), a minimal self- substitution transformation creates 
the system Eq, : y = G(")(x,y) by selecting a e [0,1] and a pair of indices i,j 
(possibly the same) with dGi{x, y) / dyj ^ and then substituting aGj{x, y) + (1 — 
a)yj for a single occurrence of yj in Gj. Suppose H{x, y^; y) is the result of replacing 
the single occurrence of yj in Gj by a new variable ayo. Then the system Eq, is 

' yi = G^r\x,y) := G,{x,y) 



G\"\x,y) := H{x,aGjix,y) + {l-a)yjy,y) 



Vm = G^m\x,y) := Gm{x,y) 

More generally, a system E* : y = G*(x, y) is a self- substitution transform of E : y = 
G(x, y) if there is a sequence Eq, Ei, . . . , E^ of systems such that E = Eq, E* = E^, 
and for < i < r the system Ej+i is a minimal self-substitution transform of Ej. 

Lemma 19. For E^ and E* as described above: 

(a) S = Eq. 

(b) // E is irreducible and a e [0, 1) then E^ is irreducible. 

(c) Suppose E is irreducible. Then E* is irreducible iff each step Ej is irreducible. 

(d) Suppose E is well-conditioned and a e [0, 1]. Then E^ is well- conditioned iff 
it is irreducible. In particular E^^ is well-conditioned if a E [0, 1). 

(e) Suppose E is well-conditioned. Then E* is well-conditioned iff it is irreducible 
iff each step Ej is irreducible. 

Proof, (a) is obvious. For (b) note that a G [0,1) implies that Gi{x,y) is still 
dependent on yj after the using the substitution yj <— aGj{x,y) + (1 — a)yj on a 
single occurrence of yj in Gi{x,y). (c) follows from noting that if any step Ej is 
reducible, then so is every step after it. (d) follows from the definition of irreducible 
and (b). (e) follows from (d). □ 

Lemma 20. Suppose S* : y = G*(x, y) is a self- substitution transform of a well- 
conditioned E : y = G(a;, y). Then the following hold: 

(a) G(a;, y) and G*(a;, y) have the same positive domain of convergence. 

(b) E* and S have the same positive solutions and the same characteristic points. 
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(c) IfTi* is well- conditioned then E and E* have the same standard solution T(x) 
and extreme point {p, r) . 

(d) //" E* is well- conditioned then the Jacobians Jg(2^, y) and JG*{x,y) have all 
entries finite at the same positive points (a,b) in the domain of G. 

Proof. It suffices to prove this for the case that S* = Sq,, a minimal self-substitution 
transform of E as described above, namely substituting aGj{x, y) + (1 — a)yj for a 
single occurrence of i/j in Gi{x, y). Let 

H{x, yo; y) = A{x, y)yo + B{x, y), 

where A{x, y) and B{x, y) are power series with non-negative coefficients, and neither 
is 0, be such that 

Gi{x,y) = A{x,y)yj + B{x,y) 
Gf)(x,y) = A{x,y){aGj{x,y) + {l-a)yj)+B{x,y). 

For item (a), ffist suppose that (a, b) G Dom+(G). Then A{a,h) and B{a,h) are 
finite, so G\°'\a, b) is finite. This suffices to show (a, b) e Dom+(G'^"^) since the other 
G^"''(x, y) are the same as those in E. Conversely, suppose (a, b) e Dom"'"(G(Q;)). 
Again A{a,h) and B(a,h) are finite, so G'i(a,b) is finite; and as before, the other 
Gj(a,b) are finite. Thus (a, b) e Donn+(G). 

For item (b) , if i 7^ j then clearly the two systems have the same positive solutions 
since yj — Gj{x, y) is in both systems. 

II i — j first note that every positive solution of E is also a solution of E^. We 
have 

Gi{x,y) = A{x,y)yi + B{x,y) 

Gf{x,y) = A{x,y){aGi{x,y) + {l-a)yi) +B{x,y) 

= A{x, y) (a{A{x, y)yi + B{x, y)) + (1 - a)y^ + B{x, y) 

= aA{x, y)\ + aA{x, y)B{x, y) + (1 - a)A{x, y)yi + B{x, y). 

Let (a, b) be a positive solution of E^. Then (a, b) solves all equations yj — 
Gj{x, y) of E where j ^ i since these equations are also in E^. Now 

h = G^r\x,y) 

= aA{a, hfhi + aA{a, h)B{a, b) + (1 - a)A{a, h)bi + B{a, b). 
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SO 



(l - aA{a,hf - {1 - a)A{a,h)Yi = (l + aA{a,h)^B{a,h). 
Since 1 + aA{a, b) is positive, one can cancel to obtain 

(l-A(a,b))6i = B{a,h), 

and thus 

bi = A{a,h)bi + B{a,h), 

which says that (a, b) satisfies the ith equation of E, and thus all the equations of 
E. Consequently E and E^ have the same positive solutions (a,b). 

To show both systems have the same characteristic points, the partial derivatives 
dGi/di/k and dCf'^ /di/k, 1 < k <m, are compared: 

dG,{x,y) dA{x,y) dB{x,y) 
dyk dyk dyk 

dGf{x,y) dA{x,y) 



dy, - dy, {»G,ix,y) + {l-a)y,) 

, ^ f dGj{x,y) \x \ ^ dB{x,y) 



so 



dGl"(x,y) aG.(x,y) _ 8A(x,y) . , 



At a positive solution (a, b) to S (hence to S*), this gives 

SO 

aG'f^(x,y) ^gi(a,b) / gG'^(a,b) 

dyk dyk V dyk 

Thus, since (a, b) is positive, one obtains jQ,(a, b) := / — Jg(«)('^) t>) from J(a, b) : = 
/ — JG{a, b) by adding a non- negative multiple of the jth row to the ith row when 
i j-i and by multiplying the ith row by a positive constant when i = j. It follows 



34 JASON BELL, STANLEY BURRIS, AND KAREN YEATS 

that det( J(a, b)) = iff det{Ja{a, b)) = 0. Combining this with the fact that S and 
have the same positive solutions shows that CV = CVa, that is, they have the 
same characteristic points. 

For the next claim, item (c), note that the composition of minimal self-transforms 
using a G [0, 1) at each step preserves the well-conditioned property by Lemma [T9l 

For a well-conditioned system E, the standard solution is the unique sequence T(x) 
of non-negative power series with T(0) = that solve the system. The standard 
solution of E is clearly a solution of Sq,. Thus if Sq, is well-conditioned then it has 
the same standard solution, and hence the same extreme point, as S, so (d) holds. 

For the final item, let (a, b) be a point in Dom"'"(G), hence a point in Dom"'"(G^"^). 
A{a, b) is finite by looking at the expression above for Gj(x, y). Then, since G^"^ (x, y) 

Gj(x, y) for j ^ i, f l33l) shows that — ^r^— — ^ is finite iff — ^-^—^ — ^ is finite, so one 

dyk dvk 

has item (e). 

□ 

Lemma 21. A well-conditioned system S : y = G(a;, y) can be transformed by a self- 
substitution into a well- conditioned system S* : y = G*(a;, y) such that the Jacobian 
matrix JG*(3;,y) has all entries non-zero. Indeed, given any n > 0, one can find a 
E* such that all nth partials of the G* with respect to the yj are non-zero. 

Proof. The goal is to show that there is a sequence Sq, . . . , of minimal self- 
substitution transforms that go from S to the desired S*, and such that each system 
Sj is well-conditioned. The following four cases give the key steps in the proof. 

CASE I: Suppose some Gi is such that all nth partials are non-zero. If Gj is dependent 
on yi (there is at least one such j) then substituting {l/2)Gi + (l/2)?/j for some 
occurrence of yi in Gj gives a well-conditioned system S' such that for G'^ = Gi and 
G'j, all nth partials are non-zero. Continuing in this fashion one eventually has the 
desired system S*. 

CASE II: Suppose — 4^ for some i. This means y,™" divides some monomial 

of Gi. Use the fact that for any j ^ i there is a dependency path from yi to yj to 
convert, via self-substitutions that preserve the well-conditioned property, a product 
of n of the yi in this monomial into a power series which has dividing one of its 
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monomials. By doing this for each j ^ i one obtains a well-conditioned G[ with 

QmuQ, 



dyi" ■■■dy, 



n 

m 



+ 0. 



Yl is now in Case I. 



CASE III: Suppose — ^ 7^ for some i. Substituting Gi for a suitable occurrence 



of yi in Gi gives a well-conditioned E' where — — ^ 7^ 0. Continuing in this fashion 

oyf 

leads to Case II. 

CASE IV: Suppose 7- — — ^ ^ for some i,j,k. U j ^ i there is a dependency 

oyjoyk 

path from yj to yi which shows how to make self-substitutions (that preserve the 

d'^G- 

well-conditioned property) leading to 7: — — ^ ^ 0. Likewise, \i k ^ i there is a 

oyioyk 

dependency path from y^ to yi which shows how to make self-substitutions (with 

each minimal step being well-conditioned) leading to -- — ^ ^ 0, which is Case III. 

oyi^ 

Since E is non-linear in y, for some i, j, k we have 

dyidyu 

Thus starting with Case IV and working back to Case I we arrive at the desired E*. 

□ 

Lemma 22. Let E : y = G(a;, y) he a well- conditioned system and let E* : y = 
G*(a;, y) he a self-substitution transform of 12. // (a, b) is a characteristic point of 
E, hence ofj:*, then A(a, b) = 1 ijf A*{a,h) = 1. 

Proof. Let (a, b) be a characteristic point of E. It suffices to consider the case where 
E* is obtained from E by a minimal self-substitution. Let Gi{x, y) depend on yj, and 
let H{x,yQ\y) be the result of replacing a single occurrence of yj in Gi{x^y) by y^. 
Then let E^ : y = G'^^\x,y), a G [0, 1], be the minimal self-substitution transform 
of E obtained by applying the substitution y^ Q;Gj(a;, y) + (1 — a)yj to H{x, yo; y) 
to obtain 

G'f'\x,y) = H{x,aGj{x,y) + {1 - a)yj);y). 
Let Aq := Aa{a, b), the largest real eigenvalue of jQ(a)(a, b). 
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The only information that we need from the above construction of the G,-"'' is that 
the function a ^— JQ{a){a, b) is continuous on [0, 1], and each JQ,{a){a, b) has 1 being 
an eigenvalue. Since A is continuous on non-negative matrices by Corollary [3l it 
follows that q; ^— s> Aq, is continuous on [0, 1]. The goal is to show that one has Aq = 1 
iff A„ = 1. 

Since (a, b) is a characteristic point of Sq it is also a characteristic point of S^, 
by Lemma [20l for a G [0,1]. Thus 1 is an eigenvalue of jQ,(o,){a,h) for a G [0,1]. 
Suppose Aq = 1. Suppose there is a /5 G (0, 1] with A^ > 1. From the continuity of 
Aa there is a 7 G [0, (3) such that: A^ = 1, and A^ > 1 for a G (7, /?]. 

Let Pa{x) be the characteristic polynomial of jQ(Q)(a,b). From 

p«(l) = p„(A„) = 

one has, for each a G (7, a G (1, A^) such that 

dpa 
dx 

Since A^ is continuous on [0, 1], 



■(c„) = 0. 



lim Aq, = A^ = 1. 

a— >7+ 

This implies 

lim Ca = 1, 

and thus 

^(1) = lim ^(c„) = 0. 
dx Q^7+ dx 

But from the Perron- Frobenius theory (see Proposition [2]) we know that A^ = 1 
implies that 1 is a simple root of p-^ix), giving a contradiction. Thus Aq = 1 implies 
A„ = l. 

A similar proof gives the converse, that if Aq, = 1 then Ag = 1, proving the lemma. 

□ 

Remark 23. In view of the last two lemmas, given a well- conditioned system E : 
y = G(a;, y), when one wants to prove something about the positive solutions, the 
characteristic points, or whether or not A(a,b) = 1 at a characteristic point (a,b), 
one can, given any n > Q, assume without loss of generality that all nth partials 
of each Gi with respect to the yj are non-zero. In the rather scant literature on 
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nonlinear systems one finds a preference for working with aperiodic systems (see, 
e.g., [7\), no doubt because of the simplicity of using uniform substitutions to convert 
such a system into one where the Jacobian of G has non-zero entries. With Lemmas 
l21\ and\2^ the need for the aperiodic hypothesis is avoided. 

4.2. Basic Properties of (p, r) and CV. Now we turn to the question of how 
to find information about the extreme point (p, t) of a well-conditioned system S 
without solving the system for the standard solution T(x). 

Lemma 24. Let y = G(a;, y) be a well- conditioned system with all entries of Jg 
non-zero. 

(a) One has the formal equality 

(34) T'(x) = G^{x,T{x)) + Jg{x,T{x)) -Tix), 

which also holds for x G [0, oo] . 

(b) All Tl{p) are finite or all Tl{p) = oo. 

(c) For all i,j the following hold: 

< 7^ P,T < oo 

< ^(p..) < 1. 

Proof. Differentiating gives ([MD, so T'(x) is a solution to the irreducible system 
u = Gx{x, T(x)) + Jg{x, T(x)) ■ u, implying (b). For x G (0,p), for each i,j, (l34l) 
implies 

T:{x) > ^{x,T{x))-T;ix), 

and thus 

1 > g(x,T(x)).^(x,T(x)) > 0, 

dG dG 
giving the inequalities in (c) since the value of -——^ (p, r) is the limit of ——^ (x, T(x)) 

as X approaches p from below. □ 
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Lemma 25. Let y = G(x, y) be a well-conditioned system. 

(a) // (a, b) G CV then A(a, b) > 1. 

(b) < A(a, T(a)) <l,forO<a<p. 

Proof. For (a) note that (a, b) G CP implies that 1 is an eigenvalue of Jg('2, b), so 
A(a,b) > 1. 

(b) Given < a < p, by the Perron-Frobenius theory of nonnegative matrices we 
know that there is a positive left eigenvector (a row vector) v belonging to A (a, T(a)) . 



Lemma 26. Let y = G(x, y) be a well- conditioned system. Suppose (a, b) and (c, d) 
are characteristic points and (a, b) < (c, d). Then (a, b) = (c, d). Thus the set of 
characteristic points of the system form an antichain under the partial ordering < . 

Proof. For the proof assume, w.l.o.g., that all second partials of the Gj with respect 
to the yj do not vanish. If b = d then G(a,b) = b = d = G(c, d), which forces 
a = c by the monotonicity of each G,. 

Now assume b 7^ d. Since b < d, all entries of d — b are non-negative. Using part 
of a Taylor series expansion, 



By O 




V ■ G,{a, T(a)) + v • J^ia, T(a)) ■ T'(a), 



so 



v-T'(a) = v-G4a,T(a)) + A(a,T(a))v-T'(a). 
Since v ■ T'(a) > and v ■ G^(a, T(a)) > it follows that A (a, T(a)) < 1. 



□ 




a^Gi{a,h) 



G(c,d) > G(a,b) + JG(a,b)(d 



a^G,n(a,b) 




Since G(a, b) = b and G(c, d) = d. 



a^Gi(a,b) 



{di - b,y 



d-b > JG(a,b)(d-b) + - 
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Let A be the largest real eigenvalue of the positive matrix JG(a,b), and let v be a 
positive left eigenvector belonging to A. Then 



v(d-b) > vJG(a,b)(d-b) + ^v 



Av(d - b) + -V 



a2G„(a,b) / , _ , 



a^Gi(a,b) 



" {dm 



SO 



:i-A)v(d-b) > -V 



dym^ 



a^Gi(a,b) 



dy 



Qy^2 'Jm) 



> 0, 



and this forces A < 1, contradicting Lemma [25] (a). 



□ 



Lemma 27. Let y = G(x,y) be a well-conditioned system. 

(a) (p, t) is in the domain of jQ,{x,y), that is, all entries of the matrix Jg{p,t) 
are finite. 

(b) // (p,t) is in the interior of the domain o/G(x,y) then it is a characteristic 
point. 

(c) < A(p,t) < 1. 

(d) A(p,t) = 1 iff 1 is an eigenvalue of Jg{p,t) iff {p,t) G CV. 

Proof. For item (a), first let S* be a well- conditioned self-substitution transform of 
S with all entries in JG*(x,y) non-zero. By Lemma [2^ all entries of Jg*{p,t) are 
finite. Then Lemma [20] (e) shows that all entries of Jg(p, t) are finite. 

For the remainder of the proof we can assume that all entries in Jg are non-zero. 
For part (b) one argues just as in the case of a single equation — if (p, r) is an interior 
point but not a characteristic point then by the implicit function theorem there would 
be an analytic continuation of T(x) at p, which is impossible. 
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For (c), since A is a continuous nondecreasing function by Corollary [3], and since 
the limit of Jg,(^x,T{x)) as x approaches p from below is Jg{p,t), it follows from 
Lemma [25] (b) that A(p,t) < 1. 

For (d), clearly A(p,r) = 1 implies 1 is an eigenvalue of Jg{p,t), and this in turn 
implies that (p, r) G CV. Now suppose that (p,t) G CV. Then 1 is an eigenvalue of 
Jg{p,t), so A(p,t) > 1. Thus (c) gives A(p,t) = 1. 

□ 

Lemma 28. Let y = G(x, y) be a well- conditioned system. If (a, b) is a character- 
istic point and (a, b) 7^ (p, r) then either 

(a) bi > Ti for all i, or 

(b) a < p and b^ > Ti{a) for all i, and some bj > Tj. 

Proof. Condition (e) in the definition of well-conditioned ensures that each Gi{x,y) 
depends on x. W.l.o.g. assume that all second partials of each Gi{x, y) with respect 
to the Uj are non-zero. Suppose that (a) does not hold. 

Claim 1: If some bi > Ti and some bj < tj then a < p and Tj(a) < bi for 1 < i < m. 
WLOG assume that 

bi<Ti,...,bk<n 

and 

bk+l > Tfe+l, . . ■ ,bm > Tm- 

From the monotonicity and continuity of the Tj on [0, p] it follows that for \ <i <k 
there exist unique ^i G (0, p] such that 

WLOG assume that 

< ei < ■ ■ ■ < < p- 

For z G {1, . . . , k} 

mi) < m^) = h 

and for k + 1 < i < m 

mi) < T,{p) <k. 

Now suppose ^1 < a. Then 

bi = mi) 
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= Gi(ei,Ti(6),...,T^(6)) 
< Gi{a,bi, . . . ,bm) = h, 

a contradiction. Thus 

< a < 6 < • • • < a < P- 
Using this one has, for 1 < i < /c: 

mi) = bi 

= Gi{a,bi,...,bm) 

fc+l? • • • 5 bm) 

> Gi{a,Ti{a),...,Tk{a),Tk+i{a),...,Tm{a)) = Ti{a). 
Thus for 1 < i < /c, 

< a < < p 
Ti{a) < Tii^i) = 6,. 
Furthermore, ior k + 1 < i < m, 

Ti{a) < Ti{p) < bi. 
Thus, in this case, for 1 < i < m one has Tj(a) < 6j. 
Claim 2: li bi < Ti{p) for all i then a < p and 6j = 7i(a) for all i. 

Choose G (0,p] such that 6^ = T^(Ci)- WLOG one can assume < < • • • < 
Cm < P- If < ^ then 

bi = Gi(a, 6i, . . . , 6^) 

= Gi{a,Ti{Ci),...,TmiCrn)) 

> Gi(ei,Ti(ei),...,T^(6)) 

a contradiction. Thus a < < • • • < < P- 
Next one has 

Cm ('Cm ! -^1 ('Cm ) i • • • i -^m ('Cm ) ) 

> G^(a,Ti(^i), . . . ,T^(C^)) 
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= Gmio.-, bl, . . . , bm) 
bm 1 

SO the > step must be an equality, and this imphes = o,. Thus all = a, and 
then for all i one has bi = Ti{a). Since (a, b) = (a, T(a)) is assumed to be a different 
characteristic point from {p,t), it follows that a < p. 

Claim 3: It is not the case that bi < ti for all i. 

Otherwise by Claim 2 we would have (a, b) = (a, T(a)) with < a < p, and then 
by Lemma 1251 it would follow that (a, b) ^ CV. But by assumption, (a, b) G CV. 

□ 

Proposition 29. A well- conditioned 1-equation system y = G{x,y) has a most one 
characteristic point; if there is such a point it must be the extreme point (p, r) of the 
standard solution T{x). 

Proof. The characteristic system is 

y = G{x,y) 
1 = Gy{x,y). 

Suppose (a, b) E CV is different from (p, r). Then 6 > r by Lemma [281 

CASE 1: Suppose a > p. Then (p, r) is in the interior of Dom+(G'), so (p, r) G CV 
by Lemma [2W b). But this violates the antichain condition of Lemma [261 for CV. 

CASE 2: Suppose a < p. Then b = G{a,b) and T{a) = G{a,T{a)) leads to 
1 = Gy{a,C,) for some T(a) < ,^ < 6. But Gy{a,b) = 1 since (a, 6) G CV, so 
again we have a contradiction by the strict monotonicity of Gy{x, y) in Dom+(G). 

Thus the only possible (a, b) G CV is (p, r). □ 

Remark 30. Meir and Moon [13j prove that well- conditioned 1-equation systems 
have at most one characteristic point in the interior of Dom^(G); and if such a 
point exists then it must be {p,t). See also Flajolet and Sedgewick [7], Chapter VII 
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Proposition 31. Suppose (p, t) is a characteristic point of a well- conditioned system 
y = G(x, y). Then: 

(a) p is the largest first coordinate of any characteristic point, that is 

p = max |a : (a, b) G Cpj, 

(b) (p, t) is the only characteristic point whose first coordinate is p. 

Proof. Use Lemmas 126) and [281 □ 

5. ElGENPOINTS 

The results developed so far do not give a practical way of locating (p, t) for 
well-conditioned systems with more than one equation. Even if one is successful in 
finding all the characteristic points, no means has yet been formulated to determine 
if (p, r) is among them. In this section special characteristic points called eigenpoints 
are shown to provide the correct analog of characteristic points when moving from 
1-equation systems to multi-equation systems. 

Lemma 32. Suppose (a, b) is a characteristic point of the well- conditioned system 
y = G(a;,y). Then A(a,b) = 1 iff {a,h) = {p,t). 

Proof. We can assume that no partial dGi/dyj is zero. The direction (<^=) follows 
from Lemma [271 (d). To prove the direction (^) assume (a, b) ^ {p,t). By Lemma 
[28] one has two cases to consider: 
(I) a > p and for all i, hi > Ti 
(II) a < p and for all i, bi > Ti{a). 

For (I), (p, r) is in the interior of the domain of G, so by Lemma [27] (b) it is a char- 
acteristic point. However this contradicts Lemma [20] which says the characteristic 
points form an antichain. 

For (II), from the equations 

G(a,b)-b = 
G(a,T(a)) - T(a) = 
one can apply a multivariate version of the mean value theorem to derive: 
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with Yij = {vij{l), . . . , Vij{m)) satisfying 

Vij{r)=Tj{a) ifr>j 
Ti{a) < Vij{r) <bi if r = j 
Vij{r) = bj if r < j. 

'da 

du 



Clearly (1351) shows that A = 1 is an eigenvalue of ( -(a, Vjj) J, and from the 



properties of the v^- we see that for all i,j 

dGi dGi 
^(a,v.,) < ^(a,b) 

since each dGi/dyj depends on all the variables x,yi, . . . , ym- 
From these remarks and the monotonicity of A one has 

1 < A(^(a,v,,)^ < A(a,b), 

showing that (a, b) ^ (p,t) implies A(a,b) > 1. 

□ 

Definition 33. A characteristic point (a, b) is an eigenpoint if A(a, b) = 1. 

The following theorem summarizes the key results proved in this section for well- 
conditioned systems. 

Theorem 34. Let S : y = G(x, y) he a well- conditioned system. Then the following 
hold: 

(a) (p,r) G Dom+(G) 

(b) If {p,t) is in the interior ci/Dom"^(G) then it is an eigenpoint. 

(c) The system S has at most one eigenpoint. 

(d) If there is an eigenpoint ofT, then it must he {p,t). 

(e) If there is no eigenpoint ofT, then {p,t) lies on the boundary o/Dom+(G) and 
one has A(p,T) < 1. 

This result can be superior to Proposition [31] for computing purposes since the 
latter requires that one know all characteristic points of S before being able to isolate 
the one candidate for (p,t). Theorem [5^ says that if one can find a characteristic 
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point (a, b) with JG(a,b) having largest positive eigenvalue 1, it is {p,t). As with 
the 1-equation case, if there are no eigenpoints of S, then new methods are needed. 

Example 35. The well- conditioned polynomial system 

yi = Giix,yi,y2) := x{l + 2y^ + 2xSly2) 
y2 = G2ix,yi,y2) := x{l + x^y2 + 2yfyl) 

has four characteristic points which, to 6 places of accuracy are: 

(0.1818598, 1.556545, 0.3647603) 
(0.2640956, 1.210710, 0.5353688) 
(0.3867644, 0.6661246, 3.834789) 
(0.4153198, 0.6217456, 0.4743552) 

One sees that these four points form an antichain, as required by Lemma The 
extreme point {p,Ti,T2) of a polynomial system is a characteristic point. By Propo- 
sition EH it must be the last one since it has the largest x-value, assuming one has 
found all characteristic roots of this system. If one is not sure that there are only 



four characteristic points then, by Theorem\3^ it suffices to verify that the indicated 
characteristic point is an eigenpoint. 



6. Drmota's Theorem Revisited 

In 1993 Lalley [10] proved that the solutions = Tj(x) to a well-conditioned 
polynomial system y = G(x, y) would have a square-root singularity at p, and thus 
one had the familiar Polya asymptotics for the coefficients^ In 1997 [5], and again 
in 2009 [6], Drmota presented the first sweepingly general theorem concerning the 
asymptotic behavior of the coefficients of solutions of a well-conditioned system, 
namely the coefficients will again satisfy the same law that Polya found to be true 
for several classes of trees (see [16]). However, as explained in Footnote Ej the 
hypotheses that Drmota has for the characteristic points of the system seem to be 
incorrect in the first publication, and vague in the second. To prove the theorem one 
needs to be able to show that (p,t) is in the interior of the domain of G(x, y). The 

^Having a polynomial system is a very strong condition since it immediately tells you that p is 
a branch point, which leads to a Puiseux expansion; it is only a matter of determining the order of 
the branch point (which is nonetheless a nontrivial task). 
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following subsection gives a clear statement of the hypotheses needed, along with a 
slightly different proof of the key induction step for the proof. 

6.1. Drmota's Theorem. The following version is somewhat simpler than that 
presented by Drmota since there are no parameters. 

Theorem 36. Let S : y = G(x,y) be a well- conditioned system with standard 
solution T(x). Suppose S has an eigenpoint in the interior o/Dom"'"(G). Then each 
Ti{x) is the standard solution to a well-conditioned 1-equation system yi = Gi{x,yi). 
Thus each Ti{x) has a square-root singularity at p, and the familiar Poly a asymptotics 
(see, e.g., [2]j hold for the non-zero coefficients. 

Proof. One only needs to consider the case that the system has at least two equa- 
tions, and one can assume all second partials of the Gt with respect to the yj are 
non-zero. The following shows that eliminating the first equation (and yi) yields 
a well-conditioned system with one less equation which has the standard solution 
(T2(x), . . . , Tm(x)) and an eigenpoint in the interior of the domain of the system. 
By the Implicit Function Theorem one can solve the first equation 

yi = Gi(x,y) 

for yi, say 

yi = Hi{x,y2,...,ym), 
where Hi is holomorphic in a neighborhood of the origin, that is, Hi{0, 0) = and 

Hi{x,y2,...,ym) = Gi{x,Hi{x,y2,...,ym),y2,---,ym) 

in a neighborhood of the origin. 

Since the Tj(x) take small values near the origin (as they are continuous functions 
that vanish at x = 0), it follows that 

Hi{x,T2{x), . . . ,r„(x)) = Gi(^x,Hi[x,T2{x), . . . , T^(x)), T2(x), . . . ,T^(x)j 

holds in a neighborhood of the origin. Also one has 

Ti(x) = Gi{x,Ti{x),T2{x), . . . ,Tm{x)) 

holding in a neighborhood of the origin, so by the uniqueness of solutions in such a 
neighborhood, we must have 

Ti{x) = Hi[x,T2{x), . . . ,Tra{x)) 
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in a neighborhood of the origin. By Proposition [H this equation actually holds 

globally for |x| < p; in particular Hi converges at (p, T2, . . . , Tm). By Corollary [3]^a) 
dGi 

the Jacobian 1 — - — of the equation yi = GUx, y) does not vanish at (p, r). Thus, 
dyi 

by the Implicit Function Theorem, Hi is holomorphic at (p, T2, . . . , Tm)- 

Now discarding the first equation and substituting Hi{x,y2, ■ ■ ■ ,ym) for yi in 
the remaining equations gives a well-conditioned system of m — 1 equations y^ = 
G*{x,y2, ■ ■ ■ ,ym), 2 < i < m, with standard solution (x, T2(x), . . . , Tm(x)) whose 
extreme point (p, T2, . . . , Tm) is an eigenpoint, since it is a characteristic point of the 
system that is in the interior of Dom'^(G*). Thus the elimination procedure can 
continue if G* consists of more than one equation. 

□ 

The extreme point of a well-conditioned polynomial system, such as Example 
[35l is always a characteristic point, and, as Lalley [lOj proved, the coefficients of 
the solutions Tj(x) have the classical Polya form Cjp""n^^/^. Drmota [5J extended 
Lalley's result to well-conditioned power series systems with the extreme point in the 
interior of the domain of the system. A natural (and desirable) direction to consider 
for further research would be to drop the irreducible requirement. However, even in 
the polynomial case, this leads to substantial challenges. For example, consider the 
reducible polynomial system 

yi = z/3 ■ (1 + 1/2 + yf) 
y2 = z/3 ■ (1 + 1/1 + 1/2) 

^3 = x-{l + 9ys^). 

Let the third equation have the standard solution y^ = A{x). One then sees that this 
example is really just an alternate presentation of Example [16] where the solutions 
for yi and y2 have a fourth- root singularity at their radius of convergence. 

6.2. A Wealth of Examples. In [2| we showed that single equation systems formed 
from a wide array of standard operators like Multiset, Cycle and Sequence led to 
square-root singularities and Polya asymptotics for the coefficients. The arguments 
used there easily carry over to the setting of systems of equations since the conditions 
in that paper force the positive domain to be an open set, and this guarantees that 
(p, r) is an interior point of the domain of the system. Thus, for example, consider 
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the collection of unlabeled rooted two-colored trees, say by the colors red and blue, 
defined by requiring that immediately below each red node that is not a leaf there 
must be a prime number of blue nodes and an even number of red nodes, and 
immediately below each blue node that is not a leaf all nodes have the same color. 
Then, by Theorem [36l the number of trees in this class has the Polya asymptotics 

7. Some Open Problems about Characteristic Points of 
Well-Conditioned Systems 

Question 1. How can one locate {p,t) if it is not a characteristic point? 

Question 2. Is the set of characteristic points always finite? 

As one has seen in the introductory examples, a system can have multiple char- 
acteristic points; the two equation polynomial system in Example [35] has four char- 
acteristic points. The following example shows that the set of real solutions to the 
characteristic system need not be finite. However Question [2] asks if the set of positive 
solutions is finite. 

Example 37. For the characteristic system (belonging to a 2-equation system) 

Vi-x- {l + yi + ym) = 

y2-x- {l+y2 + ym) = 

(x — 1) ■ (x + xyi + xy2 — l) =0 

the real solutions include the infinite curve 

{(a;, 2/1, 2/2) : X = l,2/i?/2 = -l}- 
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